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We present a novel approximation scheme for the treatment of strongly correlated electrons in 
arbitrary crystal lattices. The approach extends the well-known dynamical mean field theory to 
include nonlocal two-site correlations of arbitrary spatial extent. We extract the nonlocal correlation 
functions from two-impurity Anderson models where the impurity-impurity distance defines the 
spatial extent of the correlations included. Translational invariance is fully respected by our approach 
since correlation functions of any two-impurity cluster are periodically embedded to fc-space via a 
Fourier transform. As a first application, we study the two-dimensional Hubbard model on a simple- 
cubic lattice. We demonstrate how pseudogap formation in the many-body resonance at the Fermi 
level results from the inclusion of nonlocal correlations. A comparison of the spectral function with 
the dynamical-cluster approximation shows qualitative agreement of high- as well as low-energy 
features. 
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I. INTRODUCTION 

Compounds with strongly correlated electrons ""^'^ are 
in the focus of modern solid-state research. They usu- 
ally contain transition metal ions where part of the va- 
lence electrons move in narrow bands or are nearly lo- 
calized. Prominent examples include high-temperature 
cuprate superconductors,"^ frustrated magnets^ and 
heavy-fermion systems.^ 

Despite a large effort, various fundamental questions 
are still open.^'^"'' The model system of a single transi- 
tion metal ion immersed into a metallic matrix has essen- 
tially been solved by analytical®'^ and powerful numerical 
methods. ^"'^^ However, the more general case of regu- 
lar crystal lattices with strong electron correlations in or 
near a metallic regime is much more difficult to handle. 
The large number of relevant electronic states and the 
interplay of hybridization and interaction effects makes 
it even difficult to identify relevant degrees of freedom or 
arrive at low-energy effective models. 

An understanding of electronic correlations lies at the 
heart of most open problems in this field. Although the 
metallic character of many materials lends hope to the 
hypothesis that relevant interactions are screened and 
thus of short range, they can in principle drive corre- 
lations over all ranges of spatial distances. 

Modern approaches'^ focused on the dominant role 
of the local Coulomb repulsion. In particular, dy- 
namical mean-field theory-'^'^ (DMFT) has proven to be 
very successful. It is capable of describing lattice ver- 
sions of the Kondo effect in a regime of sizable lo- 
cal magnetic moments and provides valuable insights to 
the correlation driven Mott-Hubbard metal-to-insulator 
transition.-'^"' These approximations have met with con- 
siderable success concerning, for example, one-particle 



properties'^''^ or susceptibilities.'^ 

However, many phenomena require the inclusion of 
nonlocal correlations. Prominent examples are sys- 
tems near a quantum critical point'* or cuprate high- 
temperature superconductors."^ In recent years extensions 
of DMFT have been put forward to remedy these short- 
comings by the inclusion of some spatial correlations be- 
tween electrons.'^ In the cluster approaches, the problem 
is mapped to an effective cluster of few lattice sites and, 
in analogy to DMFT, this cluster is treated like a com- 
plex impurity in a dynamic external field. 

These approaches capture short-range correlations 
quite accurately and have contributed considerably to the 
understanding of pseudogaps and shadow bands, ^'^°"^^ 
as well as to the possible occurrence of quantum critical 
points. ^"^'^^ 

However, some problems remain in these theories. In 
cluster DMFT^s (CDMFT), for example, the transla- 
tional invariance of the crystal lattice is not fully re- 
spected. The dynamical cluster approximation^^ (DCA) 
remedies this shortcoming, but the question remains 
which choice for the size and geometry of the cluster is 
advantageous.^^ But the most fundamental limitation is 
the restriction to rather short-ranged correlations. 

In this paper we present a novel kind of self-consistent 
approach to correlated lattice systems which is in prin- 
ciple capable of including nonlocal correlations of arbi- 
trary distance. It extends the well-known DMFT by 
the inclusion of two-site correlations of all length-scales 
and is thus termed nonlocal DMFT (NLD). These cor- 
relations are incorporated by a mapping of the lattice 
model onto a multitude of two-impurity Anderson mod- 
els (TIAM) , where the impurity- impurity distance is var- 
ied. The translational invariance and crystal symmetries 
are fully respected by the scheme proposed here. 
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The details of the mapping between lattice and two- 
impurity models is presented in Sec. II. As a first ap- 
plication we study in Sec. Ill the two-dimensional Hub- 
bard model on a square lattice, a model which is usually 
considered in connection with cuprate superconductors. 
We utilize the recently developed two-impurity enhanced 
noncrossing approximation^® as the two-impurity solver 
and analyze the results obtained with respect to the qual- 
ity of our novel lattice approach and to their physical 
implications. A short conclusion and outlook is given in 
Sec. IV. 



II. SELF-CONSISTENT LATTICE THEORY 

In order to introduce the ideas underlying our self- 
consistent scheme we restrict the following discussion to 
the example of a single-band Hubbard model 

^ = E + ^fe) 4.-%.- + ^ E (1) 

k,a- j 

We use a mixed representation, where the kinetic energy 
is expressed in terms of band states with momentum k 
and spin a = {+,^} = {tii} which are created (anni- 
hilated) by the operators ^ (cj. ^r)- ^ denotes the local 
single-particle energy and the dispersion tk is the Fourier 
transformation of the hopping matrix elements Uj be- 
tween lattice sites i and j with lattice vectors and 
Rj, respectively. The lattice structure, i.e. the type of 
lattice and topology of the hopping matrix elements, is 
completely encoded into the dispersion relation t^. The 
interaction energy characterized by the Coulomb matrix 
element U is conveniently written in terms of local occu- 
pation number operators g. = Cj ^c^ ^ which measure 
an electron on site j with spin a. The lattice constant 
a = 1 will serve as the fundamental length scale and the 
spin index a will be dropped whenever possible without 
ambiguities. 

The extension of the novel approach presented here to 
more general systems with multiple bands and/or more 
complicated interactions can be readily obtained. 

The fundamental quantity of interest is the one- 
particle Green function which obeys the Dyson equation 

Gk{z) = gk{z) + gk{z)tk{z)Gk{z) (2) 

which is represented graphically in Fig. 1. The formal 
solution is given by 

Gu{z)=[gu{z)-^ -tu\'' . (3) 

The correlated part 

gk{z) = [z-e-Y.u{z)]-' (4) 

incorporates the irreducible self-energy ^^(z) and ac- 
counts for the correction to the noninteracting sys- 
tem due to the interaction term proportional to U in 
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FIG. 1. Diagrammatic representation of Eq. (2). 

the Hamiltonian Eq. (1). It describes the interaction- 
induced, i.e. correlated, part of the propagation process, 
in which single-particle transfers via elementary hoppings 
are excluded. 

It is useful to re-write the Dyson equation (2) in terms 
of the lattice T-matrix Tk , 

Gk{z) = gk{z) + gk{z)Tk{z)gk{z), (5) 

which is commonly defined via 

Tk{z) = tk + tkGk{z)tk (6) 



1 - gkiz)tk 

In the last equality the explicit form of Eq. (3) is used. 

Processes contributing to the irreducible self-energy 
'Ek{z) can be grouped together according to the nature 
and degree of correlations they contain. One contribu- 
tion represents the exact solution of an isolated inter- 
acting local site, which is represented by the "atomic" 
sclf-encrgy29 T,(o\z) = u'^^^^^^^^. Ah other terms 
incorporate genuine lattice processes. Some of these are 
captured by the well-known self-energy of the DMFT- 
approximation, which furnishes a fc-indepcndent but 
dynamic contribution ^'^^■'(z). It incorporates loops 
through the lattice attached at one site at which all lo- 
cal (dynamic) correlations are fully respected. Beyond 
DMFT, correlations between two or more of these loops 
are generated by interaction events at different sites and 
thus constitute nonlocal cumulant corrections. 

Organizing these cumulant corrections according to the 
number of different lattice sites they correlate, the corre- 
lation function of Eq. (4) can be expanded as 

5^i(z)=z-6-S(°)(z)-SW(z) (8) 
-tf{z)-tf{z)-... 

It is the aim of the present investigation to identify and 
calculate the self-energy contribution Ti),'{z), which con- 
tains all nonlocal correlations between any two sites of 
the lattice. While two-site correlations are thus explic- 
itly included in this scheme, three-site and higher-order 
nonlocal cumulant corrections E^"^(z) with n > 3 are ne- 
glected. The function incorporating these correlations, 
gk{z), will be extracted from solutions of various two- 
impurity models where the impurity-impurity distance is 
varied. Then, the Dyson equation (2), respectively its 
solution Eq. (3), is used to obtain the lattice Green func- 
tion. 
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A. Mapping between the lattice and a set of 
two-impurity models 

In order to extract the function gkiz) which includes 
two-site correlations, we consider a two-impurity Ander- 
son model'^" (TIAM). The two magnetic impurities arc 
immersed with finite distance into a host with a nonin- 
teracting conduction band. The Hamiltonian of one such 
general TIAM with a fixed distance vector a is given by 

^(2imp,^) ^ Y: .,4,.^.,. + E uiff^,..nL<r' (9) 

j={l,2},a j,l,a,(j' 
+ E ^MPi'y^.a 

+ -^E(^^^"---lA. + ^-)- 

The c-operators (d-operators) describe local interacting 
electrons on the impurities at positions Rj G {i?i,i?2} 
(the nonintcracting conduction band electrons with mo- 
mentum k) with spin a. The spatial distance between 
the two impurities is fixed and given by a = Ri — i?2- 
ej is the local single-particle energy on each impurity 

site and CZ/j^^ are the Coulomb matrix elements of the 
local (j = /) and nonlocal (j 7^ I) density-density inter- 
actions. The term proportional to t^^j^ is a possible direct 
single-particle hopping inside the two-impurity cluster 
and M^/j-^ collects additional interactions as, for exam- 
ple, exchange, pair hopping or correlated hopping. The 
last term represents the hybridization between interact- 
ing electrons of the two-site cluster and nonintcracting 
band electrons. 

The scheme proposed is capable of properly treating all 
such terms in the Hamiltonian so that nonlocal Coulomb 
interaction or long-ranged single-particle hopping matrix 
elements can in principle be considered. For the present 
case of the single-band Hubbard model nonlocal Coulomb 
matrix elements are not included and the single-particle 
hopping coincides with that of the underlying lattice as 
specified in the Hamiltonian of Eq. (1). 

Since we are only interested in properties of the inter- 
acting c-electrons neither the hybridization matrix ele- 
ments nor the dispersion of the band electrons are di- 
rectly of interest. Only the effective medium for the two- 
impurity cluster is relevant which is expressed in terms 
of these parameters as 



self-consistent theory for a lattice model a correspond- 



f{2imp.a) 1^^^ 



1 



N ^ 

k 



(10) 



It represents propagation processes via the nonintcract- 
ing conduction band in form of irreducible loops start- 
ing and ending at the two-site cluster. In our truly 



ing effective medium Ta 



(2iinp,a) 



has to be determined by 



an appropriate mapping between the lattice and the two- 
impurity model. Therefore, it becomes dressed and takes 
into account repeated interactions on local shells of other 
lattice sites. 

A Dyson equation for the local two-impurity Green 
functions is set up as follows 

gr(2imp,a) _ ^(2imp,a) j-j^j^-j 
+ g(2imp,a) f^{a) _^ ^(Simp.a) \ ^(2imp,a)^ 



where we used a matrix notation for the spatial compo- 
nents 



^(2imp,a) 



'A 



(2imp,a.) 


(2imp,a.) 



^(2imp,a)^ 



A (^zinip,fij * (2imp,a) 



(12) 



The correlated part (;(2imp,a) j-cprcsents the cluster cumu- 
lant Green function and takes the interaction matrix ele- 
ments into account. Therefore, the direct ~ t^-^ and indi- 
rect hopping events via the effective medium ^ T*' 
are explicitly incorporated in Eq (11). 

The Dyson equation (11) can be formally solved to 
yield 



(2iinp,a,) 



(^)-l_t(il)_t(2™)(^) 



£(2imp,a)(^-) ^ _ _ _ 

'(13) 

which apart from the matrix structure and the occurrence 
of the hopping matrix i*^-) has the same form as in the 
single-impurity case.^ 

It is instructive to formulate the above equations in 
terms of the full T-matrix of the TIAM, thereby estab- 
lishing an equivalence to the lattice equations (5) to (7). 
The local Green function can be expressed as 

g(2imp,a) „ ^(2imp,a) _^ ^(2imp,a)2i(2imp,a)^^(2imp,a) ^-j^^-j 

In contrast to the irreducible medium T*" the T- 

matrix incorporates repeated visits of the two-impurity 
cluster. The T-matrix is thus built up from irreducible 

loops, accounted for by the inclusion of t^-^ -I- 
for steps inside and outside the cluster, and repeated 
dwellings inside the cluster, where each visit contributes 
a factor (^(^imp.a)^ expressed as 



^(2imp,a,) 



-^(2im.p,a.) 



_ g(2imp,a) 



(15) 



A mapping between the two-impurity and the lattice 
model is accomplished by connecting the irreducible cor- 
relation functions ga^z) of both approaches. For each 

nonzero distance vector a with solution g^™^^'-\z) of 
the two-impurity model, we identify 



lg(2imp,a)(^) 



(«^Q) 



(16) 
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where the matrix elements [cf. Eq. (12)] 



.9a(^) 



—ik 



(17) 



is the Fourier transform of the lattice correlation function 
of Eq. (4) and Va the number sites which have a distance 
0'= = vSi^i to a given site. 

The prefactor of in Eq. (16) arises due to the re- 
stricted phase-space for scattering in the two-impurity 
model compared to the lattice situation. In the two- 
impurity cluster the interaction-induced scattering al- 
ways transfers the electron from one site to the other 
at distance a. In contrast, in the lattice the scattering 
process transfers an electron only in a fraction of cases 
to one specific site with distance a = \a\. 

There is some ambiguity for the correlation function 
with zero distance a = 0. While in the lattice there is 
only one such function, goiz), there exist many such func- 
tions ^^^™P'-^ from all the effective two-impurity models, 
one for each distance a. Then the question arises, which 
of all these functions should be chosen. To this end, we 
take a constructive approach, where we start from a ref- 
erence correlation function g^^^'^{z) obtained with the 
DMFT and add to it all additional correlations of repre- 
sentative two-impurity clusters at different distances 



go{z) =.9°^F^(z) 



(18) 



(2imp,a) / \ 



,DMFT 



The momentum dependent correlation function of the 
lattice theory is now given by the inverse Fourier trans- 
form 



(19) 



which directly leads to the full lattice Green function via 
Eq. (3). 

As explained after Eq (10), the effective medium 

in the framework of our self-consistent the- 
ory has to be obtained from lattice quantities. This 
is achieved by identifying the momentum dependent T- 
matrices of both models. 



(2imp) 



(^) 



>gifcay^2imp,a)^^-)_ (20) 



The momentum dependent T-matrix Tk incorporates the 
correct translationally invariant linear combination of 
two-site propagations with a fixed distance. Inverting 
the Fourier transform and using Eq. (7) we obtain the 
T-matrix for a two-impurity model with fixed distance a 



T, 



(2imp,a) 



-ik a 



1 - gk{z)tk 



(21) 



From this the irreducible medium for the two impurity 
model can be obtained by inverting Eq. (15) and inserting 
the known expression for g^'^^'^P'^) . 

Lattice and crystal symmetries are completely re- 
spected in the present approach. Lattice-translational 
invariance is incorporated by construction. Even though 
the two-impurity models are solved for distances in real- 
space, the correlation functions g^^'^P -^ (^z), which explic- 
itly exclude single-particle hopping, can be periodized for 
this purpose with the help of a Fourier transform to yield 
the lattice correlation function gk{z). Crystal symme- 
tries are also fully respected since all the two-impurity 
models for which the distance vectors a.,j are generated 
by point-group transformations, have identical effective 
media and are therefore identical. Point group symme- 
tries can even be used to decrease the computational cost 
as only one representative two-impurity model of such 
a group of symmetry-related points needs to be solved. 
The solutions of all the others follow by symmetry oper- 
ations. 

We conclude this section by summarizing the steps of 
the calculation scheme: 

1. Chose a set of Na representative distance vec- 
tors {a.j} with i = 1, . . . ,iVa, one for each set of 
symmetry-related distance vectors. 

2. Start with an initial guess for the effective media 
rj^i P:_.) each distance vector. 



3. Solve Na different effective two-impurity models, 
one for each representative distance with the 
medium P'-i)^ ^\\e result are Na (matrix) 

Green functions G^2imp,a;) ^^^^ which the corre- 
lation functions (^(^imp,^:;) g^j.^ derived via the inver- 
sion of Eq. (13), 

(2imp,a) _ "g(2imp,a) -1 _j_ ^(a) j,i'^'^™V'<LV ~^ 
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(22) 



4. Map the correlation functions of the two-impurity 
models to their lattice counterparts for fixed dis- 
tances with Eqs. (16) and (18). Fourier-transform 
these gaX^) via Eq. (19) to obtain gk{z). 



5. Use gk{z) and Eq. (21) to get the two-impurity T- 
matrix j^^™P'-i^ for each distance vector a^. A new 

guess for the effective medium is obtained 

with the inversion of Eq. (15), 



i(a) 



^(2imp,a) 



2^(2imp,a.) 
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(2im.p,a.) 



(23) 



6. Go back to step 3 and iterate until convergence is 
reached. 

One-particle spectra of the lattice are then obtained 
from Eq. (3) in momentum space and the local function 
via a Fourier transform. 
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An earlier approach of Schiller and Ingersent also uti- 
lizes a two-impurity model to extend the DMFT to in- 
clude nonlocal correlations.'^^ However there arc various 
differences to our approach. They consider an expan- 
sion in the inverse spatial dimension ^ and consequently 
employ the two-impurity model only for one specific dis- 
tance, i.e. nearest-neighbor sites only. Additionally, they 
identify the irreducible self-energies [cf. Eq. (4)] of the 
lattice and the two-impurity model. In contrast, we 
establish the mapping between the lattice and the ef- 
fective impurity models via the correlation function ga. 
This represents a crucial difference, since ga includes all 
propagation processes correlated between the two sites 
via the interaction. Thus, it also incorporates repeated 
interaction-induced scattering between two sites which 
are excluded from the irreducible self-energy but need 
to be accounted for in the lattice Dyson equation (2). 
An important point concerns the incorporation of the 
explicit hopping between nearest neighbor sites t^-\ In 
the present approach the single-particle transfers and free 
propagations through the medium are explicitly sepa- 
rated from correlation effects in both models, cf. Eqs. (2) 
and (11). This implies for the two- impurity model the 
hopping to be incorporated into the effective medium [see 
Eqs. (11) and (15)] and ensures translational invariance. 
Additionally, any imbalance in the treatment of inter- 
and intra-cluster hopping in the two-impurity model, as 
it occurs, for instance, in the CDMFT, is removed. 



III. SELF-CONSISTENT SCHEME APPLIED TO 
THE 2d-HUBBARD MODEL 

In this section we present results for the two- 
dimensional (d — 2) Hubbard model (1) obtained with 
the approximation scheme described in the previous sec- 
tion. We focus on the metallic regime at not too low 
temperatures. Long-ranged magnetic order and super- 
conducting states are thus excluded, whereas magnetic 
correlations of finite extent and the precursor regime of 
a metal-insulator transition are included and accessible. 

For the solution of the two-impurity problems we em- 
ploy a solver based on direct perturbation theory in the 
hybridization^^''^^ which is an extension of the two-orbital 
solver'^'^ and which is described elsewhere. ^^■'^^ In con- 
trast to the usual noncrossing approximation"^^ (NCA) 
this two-impurity enhanced noncrossing approximation 
includes vertex corrections which allow for the accurate 
description of finite Coulomb repulsions. 

An important requirement for the two-impurity solver 
is, it needs to be able to treat dynamic non-diagonal ef- 



fective media Ta 



(2imp,a,) 



{z). 



We study the Hubbard model on a two-dimensional 
simple-cubic lattice with nearest-neighbor hopping only. 
The half-bandwidth D = 2dt = 4t is used as unit of 
energy and we set = c ~ h = 1. The noninteracting 
spectral function of this system is shown in the inset of 
Fig. 2(a). 
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FIG. 2. (Color online) Comparison of the spectral function 
of the two-dimensional (d — 2) half-filled Hubbard model on 
a simple cubic lattice for various temperatures T and for dif- 
ferent levels of the nonlocal approximation: (a) DMFT, (b) 
NLD(l), (c) NLD(2). The parameters in all calculations are 
U/D = 1 and ^ = ~y = — -j- AH curves are calculated with 
the two-impurity enhanced noncrossing approximation as im- 
purity solver. The inset in panel (a) shows the noninteracting 
density of states. 



In order to investigate the influence of the nonlocal 
correlations we consider various stages of the scheme 
differing in the maximum distance ||a|| to be incorpo- 
rated into the approximation. This distance is mea- 
sured with a Manhattan metric indicated by ||.|| reflect- 
ing the minimum number of elementary hoppings be- 
tween the two sites. For example, the approximation 
denoted with NLD(l) includes only two-impurity models 
of nearest-neighbor lattice sites, that is, the set of dis- 
tance vectors is given by % G {(±1,0)^, (0, ±1)^}. Ac- 
cordingly, in NLD(2) this set is augmented with the four 



6 



next-nearest neighbor distances, i.e. the set of vectors is 
0,2 G {oi, (±1, ±1)^, (±2, 0)^, (0, ±2)^} and so on. 

Figure 2 displays the local single-particle spectral func- 
tion 



pH = -ilml^Gfe(c. + zO+) 

TT iV — ' ~ 



(24) 



for three different temperatures and various stages of the 
approximation. Panel (a) shows the usual DMFT so- 
lution, while panels (b) and (c) show the NLD(l) and 
NLD(2) spectral functions, respectively. The sequence 
of decreasing temperatures reveals the development of 
the well known many-body resonance at the Fermi level 
/i = a; = in DMFT [panel (a)] indicating the formation 
of low-energy quasiparticles.^^ 

The inclusion of nonlocal correlations by utilizing 
NLD(l) leads to the formation of a pscudogap inside this 
resonance [see panels (b)]. A further inclusion of next- 
nearest neighbor correlations [panel (c)] the gap widens 
and the sidepeaks become more pronounced. We cannot 
definitely decide whether or not a complete gap forms at 
zero temperature (T — )■ 0) where the spectral function 
vanishes at the Fermi energy, as it was found in a recent 
two-site DCA calculation.^^ Too low temperatures can 
not be investigated with our impurity solver due to the 
violation of Fermi liquid properties. ^^'^^ 

Additionally, in the high-energy part of the spectrum 
the inclusion of additional neighbors seems to bring out 
more of the van-Hove singularities at = ±D of the 
original unperturbed (U = 0) spectrum [see inset of 
Fig. 2(a)]. 
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FIG. 3. (Color online) Comparison of the spectral function 
for the different levels of the NLD(i) approximation. The two- 
site correlations included range from between nearest neigh- 
bor sites in NLD(l) up to all pairs of sites reachable by three 
elementary hoppings NLD(3). Other parameters are as in 
Fig. 2. 

Inclusion of additional sites with ||a|| > 2 into our 
NLD-scheme does not produce considerable changes, at 
least for this choice of parameters. This can be observed 
in Fig. 3, where spectral functions for a fixed tempera- 
ture T = 0.05D are shown for various maximal distances 



up to ||a|| = 3, i.e. including neighbors reachable by up to 
three transfer processes ~ t. The changes become succes- 
sively smaller and the curves from NLD(2) and NLD(3) 
are already almost indistinguishable. 

Our calculations clearly show the appearance of high 
sidepeaks at the borders of the pseudogap. The existence 
of sidepeaks in the low-energy spectrum is in accord with 
numerical DCA results, where, however, large cluster 
sizes where necessary to resolve gap and sidepeaks. 

We would like to point out, that the sidepeaks in the 
low-lying quasiparticle regime bear a strong resemblance 
to what is observed in a related work^^ for the TIAM. 
In our opinion we now find a coherent version of the 
splitting-scenario described there: The direct and the in- 
duced effective hopping generate bonding and antibond- 
ing molecular-like orbitals amongst different sites. The 
strength of this strongly depends on the crystal structure 
since the relevant effective matrix elements vary and os- 
cillate with distance.^* They are also very sensitive to 
the position of the Fermi level, which sets a scale for 
such oscillations. The phase-information and coherence 
of electron propagations which is necessary for such a 
splitting is fully incorporated into our approach accord- 
ing to Eq. (21). 

These aspects also manifest themselves in the sharp- 
ening of the high-energy features around u Ki D when 
going from NLD(l) to NLD(2,3). 

In addition to the bonding-antibonding splitting, an 
effective antifcrromagnetic exchange Jt = — between 
neighboring sites is generated which leads to the suppres- 
sion of Kondo-like correlations. Both effects described 
above induce a pseudogap and can lead to a splitting 
of the coherent many-body resonance at the Fermi level. 
These cases could in principle be distinguished as the 
splitting due to molecular binding should be linear in 
while the magnetic singlet-triplet splitting Jj is propor- 
tional to . However, our choice of parameter values 
implies similar magnitudes for the hopping t and the ef- 
fective exchange Jt, and both effects should be of com- 
parable size here. 

An interesting question is, whether the excitations in 
the vicinity of the Fermi surface form a dispersive band or 
are narrowly concentrated in fc-space. We show in Fig. 4 
the fc-resolved spectral function 



p(fc,w) = --ImGfc(a; 

TT ~ 



(25) 



for k along the diagonal (1, l)-direction in the Brillouin 
zone. One still recognizes the original cosine-form of a 
tight-binding band via the positions of the peaks in the 
projection onto the k — w-plane. However, these peaks 
develop a considerable width with growing |a;|, which in- 
dicates that quasiparticle excitations arc not well defined 
away from the Fermi surface. Near the Fermi level cj = 
two distinct pairs of narrow peaks are visible at wave- 
vectors k ~ ±-|(l,l)^. Even though very slight rem- 
nants of these peaks are observable at neighboring wave 
vectors too, their height is rapidly suppressed as k moves 
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FIG. 4. (Color online) fc-resolved spectral function obtained 
with the NLD(2) for wave vectors along the (1, l)-direction, 
and for e/D = -i, U/D = 1 and T/D = 0.05. Panel (a) 
shows the full energy interval while panel (b) displays a mag- 
nification of the energy interval around the Fermi-energy. 

away from ±^(1,1)'^. Since these maxima apparently 
furnish the spectral weight of the flanks of the pseudogap 
observable, e.g., in Fig 3, we conclude those side-peaks to 
be rather localized in fc-space. Thus, they are a result of 
increasing lattice coherence at decreasing temperature. 

Remarkably, also near the band edges at fc = and 
k = ±7r(l,l)^ with respective energies lu w ztD rather 
narrow peaks appear on top of the broad resonances con- 
nected with the original band. Since they show up only 
if next-nearest neighbors are included, i.e. in NLD(2) 
but not in NLD(l) (and DMFT), we attribute them to 
next-nearest neighbor correlations, possibly of magnetic 
nature. A participation of an indirect binding effect due 
to the repeated action of the transfer t, i.e. a molecu- 
lar orbital-like effect, can also not be excluded with our 
choice of parameters. 

Up to now a symmetric situation 2e + U = was in- 
vestigated, in which a possible Mott-Hubbard gap as well 
as a possible pseudogap in the low-energy quasiparticle 
regime both were to open around the Fermi-level. Both 
effects can be separated and identified individually by 
moving away from half-filling, i.e. with increasing dop- 
ing. Equivalently, we increase the Coulomb-repulsion U 
beyond the value U = —2e and keep all other parameters 



fixed. 
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LO/D 
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LU/D 
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FIG. 5. (Color online) Spectral functions for fixed e = — -j 
and temperature T/D = 0.05 for various U. (a) DMFT (b) 
NLD(l). (c) Spectral function for fixed U/D = 1.5 and dif- 
ferent stages of the approximation as indicated. The inset 
shows a close up of the low-energy region around the Fermi 
level. The fillings corresponding to the curves shown in (c) 
are tidmft ~ 0.90 and TiNLofj) ~ 0.92. 

Calculations for such particle-hole asymmetric situa- 
tions are shown in Fig. 5. Panel (a) displays the forma- 
tion of a Mott-Hubbard gap at positive energies with in- 
creasing U as calculated with the DMFT-approximation. 
As expected, the center of the incipient gap moves away 
from the Fermi level by an amount proportional to U 
while the many-body resonance remains pinned at the 
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Fermi level. 

Figure 5(b) shows results of NLD(l)-calculations for 
the same values of U . The Mott-Hubbard gap forms at 
positive energies in a similar fashion as in the DMFT- 
rcsults of panel (a). Additionally, a pseudogap emerges 
in the many-body resonance at the Fermi level as it was 
already found in earlier woik?'^^'^^'^^ 

The last part (c) of Fig. 5 shows the effect of including 
more neighbors in our NLD(i)-scheme for fixed Coulomb 
interaction U . One again recognizes the dominant role of 
nearest-neighbor correlations since the spectra with i> 1 
all differ considerably from the DMFT-result. The low- 
energy region seems to be converged for i = 2, and the 
inclusion of next-nearest neighbors apparently smoothens 
the spikes around the pseudogap. 

In contrast, the inclusion of next- next-nearest neigh- 
bors {i — 3) still has considerable impact on the high- 
energy features at negative energies while at positive en- 
ergies only small differences are induced. This appar- 
ently indicates a reduced hybridization and correspond- 
ing band narrowing for states in the vicinity of the neg- 
ative band edge. The importance of lattice structure, 
coherence and correlations for the possible development 
of hybridization gaps and band splittings as they were 
already discussed in the symmetric case is evident from 
these results. 

The longer-ranged correlations also seem to work in 
favor of restoring particle-hole symmetry in the low- 
energy region [see inset of Fig. 5(c)], a phenomenon 
which is also discussed in connection with the cuprate 
superconductors. ^^'^^'•^^ 



The qualitative features of both approaches nicely agree: 
the humps and minima visible in the DCA spectra trans- 
late to corresponding but more pronounced features in 
the NLD(3). Given the different nature of the two ap- 
proximations, this provides additional evidence for the 
physical nature of the observed structures. 

The magnitude of the pseudogap is of the same order 
in both approaches, although the DCA result is slightly 
larger and the low-energy spectrum much more asymmet- 
ric. One reason for this can be found in the enhanced two- 
impurity noncrossing approximation as the two-impurity 
solver, which is known to produce a too small many- 
body low-energy scale for the SIAM. This is especially 
relevant for the low-energy region and also translates to 
self-consistent approximations (see, e.g. Ref. 15). 

The different temperatures of both calculations should 
not be relevant for the qualitative features as T is already 
rather low and further decreasing it will only slightly 
deepen the pseudogap. 

However, a qualitative difference between the two cal- 
culations is found in the different fillings which in gen- 
eral strongly influences the low-energy spectral function. 
The analysis in Ref. 24 revealed that decreasing the fill- 
ing from riDCA ~ 0.95 to tidca ~ 0.88 (for otherwise 
identical parameters) leads to a spectral function which 
is nearly particle-hole symmetric at low energies and does 
not exhibit a pseudogap. In this light, the NLD(3) spec- 
tral function for nNLD(3) ~ 0.92 with a nearly symmetric 
low-energy spectrum and a small pseudogap represent a 
plausible intermediate solution. 



NLD(3) n = 0.92 
DCA n = 0.95 




FIG. 6. (Color online) Comparison of NLD(3) and DCA 
spectral functions. The NLD(3) curve is the one shown in 
Fig. 5(c). The DCA result is taken from Fig. 5 of Ref. 24 
and was obtained from self-consistent quantum Monte Carlo 
solution of a 16-site cluster for U/D = 1.5 and T/D = 0.017. 
The filling for the DCA calculation of nocA ~ 0.95 is slightly 
larger than the NLD(3) filling, nNLo(3) ~ 0.92. The inset 
shows an enlargement of the region around the Fermi level. 

Figure 6 compares a NLD(3) spectral function for finite 
doping and U = 1.5Z? with a spectrum obtained from a 
DCA calculation as it was published in Fig. 5 of Ref. 24. 



IV. CONCLUSION 

We have proposed a novel general and nonperturbative 
scheme for the treatment of correlated electron systems 
on crystal lattices. Our approach represents an extension 
of the well-known DMFT to include nonlocal two-site 
correlations with arbitrary spatial extent. 

The self-consistent formulation establishes a map- 
ping between nonlocal two-site correlations of the lat- 
tice model and the equivalent functions of various two- 
impurity Anderson models with varying distance. The 
two-impurity models are solved in real-space, but the ex- 
tracted correlation functions arc transformed into mo- 
mentum space via a Fourier transform. Thereby transla- 
tional invariance is built into our scheme by construction. 

Most important, the spatial range of correlations which 
are explicitly included in the treatment is unrestricted 
in principle. It corresponds to the maximum impurity- 
impurity distance in an effective TIAM which is solved 
in our approach, and is therefore only limited by the ac- 
curacy of the two-impurity solver. 

As a first application, we applied our scheme to the 
Hubbard model on a two-dimensional simple-cubic lat- 
tice. We found that the leading nonlocal correlations pro- 
duce a pseudogap in the low-energy single-particle exci- 
tation spectrum as it is to be expected. ^^"^^''^^ Moreover, 
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pronounced side-peaks occur as a signature of increasing 
coherence. In a situation without particle-hole symmetry 
we could clearly discriminate between the Mott-Hubbard 
gap induced by mostly local atomic correlations and the 
pseudogap, which is found in the many-body resonance 
at the Fermi level. We have also demonstrated that the 
inclusion of correlations over larger distances brings out 
more details of the excitation spectra in the high-energy 
region. 

Additionally, our method compares very well to results 
from the DCA. The spectral functions of both approaches 
exhibit similar qualitative features for finite doping, al- 
though they are more pronounced in our scheme. Also, 
the pseudogap obtained within both approaches is of the 
same magnitude. 

Thus, the scheme presented opens an excellent per- 
spective for more detailed investigations of systems where 
nonlocal correlations play an important role. These in- 
clude phase transitions and critical phenomena in gen- 
eral, and quantum critical points in particular, as they 
are found, for example, in heavy fermion compounds or 



cuprate superconductors. 

Also, the calculation of nontrivial critical exponents 
and their scaling behavior^*^ seems to be in reach using 
this new approach. In some cases it will nevertheless be 
necessary to include higher-order irreducible correlations 
beyond those between only two sites. This will be of 
particular importance, when complicated ground states 
involving correlations on plaquettcs of sites or general 
resonating valence bond-states are under consideration.*-'^ 
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